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In the present work, we investigate iiow structural defects in graphene can change its transport 
properties. In particular, we show that breaking of the sublattice symmetry in a graphene monolayer 
overcomes the Klein effect, leading to confined states of massless Dirac fermions. Experimentally, 
this corresponds to chemical bonding of foreign atoms to carbon atoms, which attach themselves 
to preferential positions on one of the two sublattices. In addition, we consider the scattering off 
a tensor barrier, which describes the rotation of the honeycomb cells of a given region around an 
axis perpendicular to the graphene layer. We demonstrate that in this case the intervalley mixing 
between the Dirac points emerges, and that Klein tunneling occurs. 
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I. INTRODUCTION AND MOTIVATION 



Graphene, a monatomic layer of carbon atoms on a honeycomb lattice, has been synthesized for the first time 
in 2004^. Its remarkable properties, such as unconventional quantum Hall effeclj^i^, Klein tunneling^i^ or charge 
confinement^, are mainly a consequence of the fact that at low energies, the charge carriers in graphene are described 
by the relativistic Dirac equation instead of the more familiar Schrodinger equation^i^. The single particle dispersion 
relation is linear in the momentum, Ek ~ ivF\k\, where vp is the Fermi velocity and k ~ {kx,ky) is the fermion 
momentum, measured relative to the incquivalent corners of the Brillouin zone K and K' , known as Dirac points. The 
peculiar properties of graphene triggered immediately a lot of interest, due to its possible applications in nanoscale 
devices, but also because graphene is able to connect different branches of physics. Meanwhile, we have learned how 
to produce multilayers of graphene, have started to understand the effects due to the deviations from a perfect carbon 
honeycomb structure and, most important from the point of view of technological applications, are learning on how 
to grow samples with structural defects. Defects can be used to tailor graphene based devices — see, for example, the 
discussions below. However, they can also change the transport properties of the charge carriers and be at the origin 
of the performance deterioration of these devices. Comprehensive overviews of graphene can be found in Refs. Il0l4l4l 
and the structural defects are reviewed in^. The work reported in^ gives an example of an experimental realization 
of one dimensional (ID) defects in pure graphene. 

Typically, graphene based devices require the ability to confine and control the charge fiow. Confining Dirac 
fermions is yet a challenging task, due to the so-called Klein effecti^rJ^, where a relativistic electron is able to 
penetrate a barrier higher than the electron's energy completely unrefiected — in contrast to the conventional tunneling 
where the transmission probability drops off exponentially as the barrier gets higher. Studies of various types of 
barriers that lead to Klein tunneling for electrons in graphene can be found in the literature, and experimental 
observations of this effect have been reported^Sr— . Despite the difficulties related to overcoming the Klein effect, 
charge confinement in graphene has already been accomplished experimentally — see for example Refs. [7l ll0ll26l for 
overviews. Theoretically, confinement was demonstrated by cutting or bending graphene sheets^ii^^, exploiting the 
transversal degrees of freedom of the electrons in an electrostatic potential2£, applying magnetic fields^i^, deforming 
the graphene membranes^ or by spatial modulation of the Dirac gap^^*^. From structural defects, confinement has 
been achieved in graphene-graphane systems, see for example^. However, so far it has not been explored a way 
to produce charge confinement in association with chemical bonding of foreign atoms to carbon atoms, although 
techniques to implement the chemical bonding of some adsorbents, such as fluorine, hydrogen or oxygen, already 
exisli^^i^. 

In addition to the possible technological applications, the ID defects in pristine graphene are interesting because 
they can be modeled as potential barriers, associated with different fermionic operators in the Dirac equation. Within 
the conceptual framework of effective quantum field theories, the authors oS^ generalize the results o^^^ — where 
vortex formation in graphene is described through a chiral gauge theory — and introduce a scalar and a gauge field 
that account for the dynamics of the self-interaction of the carbon background and the mean self-interaction of the 
Dirac fermions. This model is well suited to describe various disorder phenomena such as topological defects, doping 
defects or distortions of the lattice honeycomb. Within this framework, carbon nanotubcs and graphene have been 
studied and their quantum properties have been reproduced^^ii^. 

In the present study, inspired by the work of^^— , we explore how charge confinement and Klein tunneling can 
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be induced by certain types of defects, and examine how defects — modeled as ID potential barriers — can be 
mapped into fermionic operators. As described below, charge localization can be achieved via a barrier which breaks 
the sublattice symmetry. In practice, the breaking of the sublattice symmetry can be realized by binding covalently 
foreign atoms with particular carbons, whereas for the theoretical analysis, we have to identify and explore the 
corresponding fermionic operators. The Klein effect is investigated by considering the scattering of an electron off 
a tensor barrier. The associated fermionic operator generates a rotation around the z-axis (perpendicular on the 
graphcnc plane) and it couples both the two sublatticcs A and B, and the two valleys K and K'. Experimentally this 
can be implemented by topologically distorting the graphene layer via a rotation of the honeycomb cells around the 
z-axis in a particular region, relative to the remaining graphene sheet. In order to include the effects of the K K' 
mixing, we work in the four-component, as opposed to the standard two-component formalism. We demonstrate that 
for this type of barrier a peculiar effect arises, i.e. the splitting of the electron wave function into two components 
inside the barrier, and we analyze its implications for the solution of the tunneling problem. 

The paper is organized as follows. In Sec. II we briefly review the relevant features of the gauge model proposed 
in Refs. I39ll40l . concentrating on the operators that account for the effects of the carbon background and the carbon- 
fermion interaction. In this setting, we construct different types of barriers associated with defects in graphene, and 
discuss the implementation of the corresponding operators in the Dirac equation. Sec. Ill is dedicated to charge 
confinement. The transmission probability of an electron scattered off a barrier with sublattice symmetry breaking is 
investigated, and the conditions that enable the clustering of charge are identified. In Sec. IV, we focus our attention 
on the Klein tunneling, and consider the scattering off a barrier that describes a spatial distortion of the graphene 
sheet. As before, we discuss the angular behavior of the transmission coefficient, underlining the effects related to the 
mixing of the Dirac points. In Sec. V a short summary and the conclusions will be presented. 



II. THE THEORETICAL SETUP 



Let us begin by sketching the fermionic dynamics in graphene-based materials. Tight binding models provide a 
first approximation to describe the electronic properties of graphene, as they only take into account the fermionic 
degrees of freedom, whereas the carbon interactions are summarized via the hopping parameters which control the 
electron dynamics. Since the dynamical degrees of freedom in graphene must include both the fermions and carbon 
background, a relativistic-like gauge model for graphene and nanotubes is suggested, which takes into account 
both the electron-hole and the carbon dynamics. In this framework, electrons are described by a four component 
Dirac- type spinor, while the carbon degrees of freedom are associated with a scalar field f and a gauge field A^. Since 
the full description is somewhat long-winded, here we briefly present the main features of the model and direct the 
reader to^ for a full account. In this gauge model, the electron dynamics is described by the Dirac equation 

{irD^-P{^)^P^{^)^5}^ = Q, (2.1) 

where the polynomials P{lp) and Psi^p) describe the interaction between fermions and the carbon crystal structure. 
Explicitly, the vacuum expectation value of the scalar field (p vanishes for pure graphene, while for doped graphene 
{(f) ~ (fo 0, where (po is a minimum of the scalar potential V{(p^ip). If one is able to fabricate graphene with islands 
where it is doped, one can simulated its electronic properties by taking 

{0 for undoped regions, 

(2.2) 
(52 + i /i2 75) 'yi'o for doped regions 

in the Dirac equation, where g2, /i2 are the corresponding coupling constants. In this way, one reduces the problem to 
the investigation of the solutions of the Dirac equation in a square potential barrier, which includes the contribution 
Eq. (|2.2p and can in principle account for any combination of Dirac 7 matrices. Furthermore, given that the potential 
barrier is associated with the doping of graphene, experiment can help to calibrate the numerical values of 52'y3o and 
h2'PQ- The covariant derivative is ~ df^ + igA^, and the four-component Dirac spinor is given by 

r+ 

V'- 

where the upper index refers to the sublatticcs A or B and the lower index to the Dirac points K (plus sign) and K' 
(minus sign). For the Dirac 7 matrices we will use the representation 
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Operators 


(S,S) 


/, 70, 7:1:75 , Sa: 


(S,A) 


75, Ix, (^Ox 


(A,S) 


Jy, 72, ^Oy, O-Qz 


(A,A) 


7y75, 7^75, ^y, 



TABLE I: The first column refers to tfie operator symmetry under intercfiange of sublattice and intervalley indices, respectively. 
The table assumes that the grapliene sheet is on the x — y plane. The E matrices are defined as a'-' — — ie'-'''"Efe, where e'-''° is 
the three dimensional Levi-Civita, e^^^ = 1 and a sum over the index k is implicit. 



where cr^ stands for the j Pauli matrix. Since we are concerned with two dimensional fcrmions, the index fi takes the 
values /X = 0, 1, 2. The Dirac matrices satisfy {7^, -f^} — 2g'^^. 

In order to simulate other dynamical effects, which can be associated with the carbon structure or with the carbon- 
electron interaction, or to include the effects due to impurities/defects, new operators should be added to the Dirac 
equation (|2.ip . These operators can be, in principle, any of the remaining matrices 7^, 7''75, cr^'^ (in this case the 
values of extend to the full four-dimensional Minkowski space). The symmetry properties of all the 16 operators 
under interchange of sublattices and of the Dirac points, respectively, are listed in TableU] One can proceed and provide 
an interpretation for each of the 16 operators of the Dirac algebra in terms of graphene defects and corresponding 
barriers. In the following, we will not provide the detailed form of each contribution; instead, we will concentrate 
on the operators that have a simple physical interpretation, and represent the experimental situations considered in 
this work. If one adds a new term to the Dirac equation, let us say —Tip, it contributes to the Hamiltonian of the 
system as ip Tip, which can be rewritten in terms of the spinor components V'+j ■ • ■ For the particular choice of F = 7^ 

(assuming that x ~ y is the graphene plane) it follows that ipj^^p = [("0+)^^+ + (V'-)^''/-] ~ [(V'+)^V'+ + 
This term distinguishes the two carbon sublattices, it is antisymmetric under interchange of the sublattice indices, 
and favors the occupation of sublattice B relative to sublattice A. Of course, if the coupling constant associated with 
this operator is negative, the occupation of sublattice A is preferred. As discussed previously, the breaking of the 
sublattice symmetry can be realized by binding covalently foreign atoms with particular carbons in graphene and, 
in this way, a 7^ operator can be simulated. It turns out that this operator (and its connection with ID potential 
barriers) overcomes the Klein effect and leads to charge confinement. A second operator of interest for this work is the 
tensor interaction 7^7^ = iSz, where Sz is the third component of the spin operator (perpendicular on the plane of 
the graphene). It describes the rotation of the graphene around the z axis and can be related to a topological defect 
of the lattice in graphene (see for example^ for a study of the low energy properties of graphene due to distortions 
of graphene sheets). Importantly, the aforementioned operators change the fcrmion dispersion relation, such that the 
electrons acquire an effective mass, with an energy gap twice the effective mass. For an explicit derivation of the 
fermion gap, and its correlation with the underlying scalar field, we refer the reader to Ref. l39l . 

Since the coupling with the scalar field ip devise the electronic properties of the doped regions of graphene, in the 
following we will set the gauge field A^ = 0. After adding the relevant operators to the Dirac equation ()2.ip . we 
obtain: 

{t^d^ -Ds^ iDpf - Doj" - DvY' - ^i?T7'7'} ^ = 0, (2.5) 

where we have introduced the notations Ds = Pif) and Dp = —iP^{ip) for the scalar and pseudoscalar interactions, 
respectively. Dj-^y^ denotes the strength of the tensor (vector) interaction, and the term proportional to 70 can be 
viewed as a chemical potential. Starting with the above equation, the task is now to investigate the transmission of 
the electrons through a ID barrier, and establish under which conditions the charge confinement and Klein tunneling 
take place. 

III. CHARGE CONFINEMENT 

In this section we restrict the analysis to a set of operators that include the scalar, pseudoscalar and vector barriers 
in the Dirac equation. Setting Dt = 0, our starting Dirac equation takes the form 

{i-f^'d^-Ds-iDp-f^ ^Do-/° -DvY}^^0- (3.1) 

The scattering over an electrostatic barrier (within the two-component formalism) was already investigated in 
Ref. for single layer and double layer graphene, and the Klein paradox has been demonstrated. When one takes into 
account in the Dirac equation the pseudoscalar interaction ipj^ip, the sublattice symmetry breaking operator tpj^^p 
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FIG. 1: Sketch of a one-dimensional defect in an infinite slab of graphene. 6 and 9' axe the angles used in the scattering problem 
in regions I (pure graphene), // (graphene with defect), /// (pure graphene), and L is the width of the barrier. 



and the scalar term '0 ijj, the scattering pattern computed in Ref. changes dramatically. In particular, the tunneling 
associated with the Klein effect can be avoided, leading to charge confinement. 

Following^, we will consider the simple geometry of an ID defect in graphene, where the doping occurs only for the 
region for < a; < i. The "experimental" setup is illustrated in Fig. [TJ The wave function describing an incident 
electron with energy E and momentum p = EicosOex + sin^Cy), in the region / (a; < 0), is given by 



where the momentum p' = E{— cos 9 
region /. For < a; < L one has 



/ Ak \ 




{ Rk \ 






Rxe-^' 






-Rk'c'" 


V Bk' ) 




\ R^' J 



(3.2) 



sinScy) corresponds to the solution associated with the reflected wave in 
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ck 
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E-Do±Di 
Ds + iDp 



and 



V„ = — 
' D 



k1 



I A-' 



V 





1 



(3.3) 



iDf 



(3.4) 



As before, q~ |(7|(cos6'' e.^ + sin6'' Cj,) and tf' = |(7|(— cos 6*' e^; + sin0' ej,). The Dirac equation has a nontrivial solution 
in region // if the corresponding determinant vanishes, which leads to the condition 



{E - D^f -Dl-Dl- D%. 



(3.5) 



As shall shortly become clear, it is important to notice that the component = \/ cp — E'^ sin^ 9 can be real or a 
pure imaginary number, depending on the parameters entering Eq. (|3.5p and the angle of incidence. For a real qx the 
Klein tunneling is allowed and the barrier is perfectly transparent for certain angles of incidence. On the other hand, 
when qx is imaginary the barrier becomes opaque and transmission is highly suppressed. For x > L the wave function 
is given by 



Tpiiiix) ^ e 



ip-x 



tK 

tK'C-' 
tK' 



J 



(3.6) 



For a relativistic system the wave function is continuous everywhere, therefore tpii^ = 0) = ipnix = 0) and 
ijjn{x = L) = ijjiii{x = L). From the matching conditions and with the definition of the transmission probability 
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FIG. 2: (Color online) Angular behavior of the transmission probability T{ff) for a 7^ barrier (blue line) plotted against a 7" 
barrier (red line), for the incoming energy _B = 80 meV and different heights of the barrier. The widths of the barriers are 
1/ = 110 nm (upper panel) and i = 50 nm (lower panel). In the figures, Vb denotes the height of the barriers, i.e. Vb = Do 
for the 7'' barrier and Vb = Dy for the 7" barrier (in general, in the second case Vb = a/ D\ -\- Dp + Dy if one takes into 
account all barriers associated with graphene doping: see also text for details). The 7° barriers reproduce the results found in 
the two-component formalism, as described in Refs. [aliol. 



one obtains the following expression: 

^ g^cos^gcos^ 9' 

^' ^cos^{q^L)cos^ecos^e' + sm\q.^L){E- Do-l^smesme')"^' 

where 9' ~ aTcta'n{py / q^) ■ In the above formula, we notice that the the Dirac points K and K' do not mix, in the 
sense that terms of the type AkB'^, are not present. Further, note that the transmission probability corresponding 
to the 7" operator, with Ds = Dp = Dy = 0, reproduces the solution obtained in Ref. within the two-component 
formalism (stemming from the four-component formalism, when the couplings between components are discarded). 
We remind the reader that in this case the result (13.81) reduces to: 



To{9) ^ ^u^™^^ 

cos2((7^L)cos2 6icos2 6'' + sin2(g^L)(l - ssin6lsin6")2 

where s = sgn{E — Dq). 

Since the 70 interaction is not relevant for our purpose to demonstrate charge confinement, in the following we will 
set Dq ~ 0. The potentials Ds, Dp and Dy all generate an electronic massive- type dispersion relation, see Eq. p.5p . 
In this sense, the doping associated with any of the scalar, pseudoscalar and 7^ operators change the inertia of the 
charge carriers in a similar way. More, given that the /, 75 and 7^ contributions in Eq. ()3.5|) are symmetric, it is 
enough to study one of the three barriers, say that associated with the 7^ interaction Dy (which, as discussed, breaks 
the sublattice symmetry), while setting Ds and Dp to zero. Then the condition Eq. p.Sp simplifies to 

f = - Dir. (3.10) 

The apparent singularities in the wave function, Eq. p.3p . can be discarded simply by multiplying the spinor with 
a factor {Ds + iDp), which changes its length but does not affect the physical results. Moreover, we have checked 
that setting Ds = Dp = Dy = directly in the Dirac equation leads to the transmission coefficient To{9), Eq. p.9p . 
although the solutions of the Dirac equation inside the barrier look rather different. 

In Fig.[2]the angular behavior of the transmission probability for a 7^ barrier, Eq. (|3.8p . with if given by Eq. (|3.10p . 
(blue line) is plotted against the 7° barrier, Eq. (|3.9p (red line), for various heights and for two barrier lengths. As 
shown in the left panels of Fig. [21 for a 7^ barrier lower than the energy of the incident electron, there are several 
directions for which the barrier is transparent. This is the "standard" Klein tunneling — already demonstrated for a 
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7o barrieriS^ — which occurs due to the fact that the wave function inside the barrier, Eq. (|3.3|) . has an oscillating 
behavior (the momentum |g| is real in this case), just like the wave function outside the barrier. Considering normal 
incidence [6 = 9' = 0) and keeping Dy < E, we notice that the 7^ barrier is not perfectly transparent, as opposed 
to the 7" barrier. Instead, transparency occurs only for higher angles of incidence. On a 7*^ barrier, for = 0' = 0, 
To{0) = 1 for any value oi QxL, independent of the height of the barrier, as can be seen from Fig.[2]and from Eq. p.9p . 
whereas in a 7^ barrier, r(0) is no longer unit but depends on the ratio E/Dy- A straightforward calculation gives 



r(o) 



1 



(3.11) 



and therefore T(0) < 1, for Dy < E. Increasing the barrier height, the range of angles where the perfect tunneling 
through a 7^ barrier is allowed shrinks until the Klein tunneling is suppressed, as can be seen in the last plots of Fig. [5] 
This behavior can be understood by looking at Eq. p.Sp . Indeed, setting E = Dy, it follows that (f = and T{0) = 
for any angle of incidence. When the height of the barrier is higher then the electron energy, the momentum q^, 
Eq. p.lOp . flips the sign such that and the component q^ become imaginary, and the wave function inside the 

barrier switches from oscillating to exponential decay. In turn, this implies that for barriers higher than the electron 
energy, the barrier becomes opaque, thus leading to the charge clustering of the electrons. 



IV. KLEIN TUNNELING AND INTERVALLEY MIXING 



So far, in the scattering off a barrier of scalar, pseudoscalar and vector type, the transitions K K' did not occur, 
in the sense that a term that mixes the coefficients corresponding to the two different valleys was not present in the 
transition amplitude. In the present section, we investigate the effects on Klein tunneling emerging from the K K' 
transitions, by introducing a tensor term proportional to 7^7^ = iSz (and discard all other contributions), where S2 
is the third component of the spin operator in the four-spinor representation. Recall that such an operator can be 
associated with a rotation of the honeycomb carbon lattice. In this case, the Dirac equation (|2.5p takes the form 



(4.1) 



As before we will consider the geometry of an almost ID defect in graphene, such that the electron is scattered off a 
barrier that ranges between and L (see also Fig.[T]). For an incident electron moving with energy E and momentum 
p — E{cos9ex + sm0ey), the wave functions in the regions / (x < 0) and /// {x > L) are given by Eqs. (|3.2p and 



p.6p . Inside the barrier, where the contribution of the tensor interaction is non- vanishing, a plane wave nontrivial 
solution for the Dirac equation requires 



D^^ = iE±\q]r. 

For a given electron energy E > there are two possible momenta, 
given by, respectively. 



|g<l,2)| 



(4.2) 

\E zt Dt\, with the associated spinors 



(4.3) 



This means that, in the scattering through the ID tensor barrier, the incoming wave is divided inside the region II, 
see Fig. [TJ into two distinct plane waves, propagating with different momenta and at different angles, 0^^^ and 0^^\ 
Since one wave is exponentially enhanced and the other, exponentially suppressed, their mixing will produce a pattern 
that exhibits perfect tunneling for several angles, regardless of the height of the barrier. This is in contrast to the 7^ 
type of defect considered in the previous section, where for a barrier higher than the electron energy the transmission 
probability goes very rapidly to zero. Using Eq. (|4.3p . we can write down the general solution of the Dirac equation 
for < a; < L: 
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(4.4) 
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where q^^^'^^ ' ' 

As in the previous section, we use the fact that a solution of the Dirac equation is a continuous function everywhere, 
in order to obtain the boundary conditions for the wave function ip. For the tensor barrier, it is convenient to write 
the boundary conditions at the borders of region II in matrix form. At .t = we have 
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With the notations 
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(4.6) 



we can rewrite the above equation as 
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V Rk' J 




(4.7) 



Similarly, the boundary condition in x ^ L reads: 



M 



tK \ 




1 e^€'^ 


tK' 







TK 











I 
























(4.8) 



Notice the insertion of the coefficients rK,rK' in the transmitted wave. Even though this alteration does not change 
the final result, there are two reasons that motivate it. Firstly, with this modification the matrix that multiplies the 
outgoing spinor becomes identical to the matrix M in Eq. (|4.7p . corresponding to the incoming wave — otherwise, 
if = rfc' = the matrix AI should be modified such that the last two columns are zero, although formally they 
can be left unchanged since the multiplication with zero does not change the final result. A second (physical) reason 
is related to the possibility of investigating a double (or multiple) square potential, where the wave emerging from 
the first barrier docs have a reflected component, which is then scattered on a second barrier, and so forth. With the 
notation 



L 
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(4.9) 



and putting Eqs. (|4.7|) . (|4.8|) together, we find that the relation between the coefficients corresponding to the incoming 
and outgoing waves can be written in matrix form as (here we set rK,TK' back to zero and remove the phase factor 
gip^L^ which does not contribute to the transition amplitude) 




= M-^NLN-^M 




(4.10) 



Furthermore, in order to calculate the transmission probability, one can rewrite the above equation in a form that 
eliminates the dependence on the reflection coefficients Rk,Rk'- With the notation 
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(4.11) 
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FIG. 3: (Color online) Angular behavior of the transmission probability T{9) for a tensor barrier 7^7^ = iEz (blue line) plotted 
against a 7" barrier (red line), for the incoming energy E — 50 meV and different heights of the barrier. The widths of the 
barrier are L — 110 nm (upper panel) and L = 50 nm (lower panel), and Dt denotes the height of the barrier. 



one can split the equation (|4.10p into "blocks" , and rewrite it in the following form: 

' ^" ( Bk' ) + 

= X2I ( n^' I + X22 



tK' 



Rk 
Rk' 



Rk 

Rk' 



(4.12) 
(4.13) 



where Xij are 2x2 matrices, corresponding to the four entries of the matrix X, Eq. (|4.1ip . It is straightforward to 
show that 



tK' ) " ^^'^ " ^12^2/^21) ( 



(4.14) 



and hence we have an explicit relation between the coefficients that contribute to the transmission probability. We 
have used Mathematica to tackle our problem and the calculation shows a matrix Xn with non- vanishing off diagonal 
terms, i.e. that produce a mixing between the coefficients Ak and Bxi. This implies that a rotation of the honeycomb 
carbon lattice induces an intervalley transition (K ^ K') in the barrier reflection and transmission. The various 
momenta and angles entering the above equations are given by 
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: EcosO, py — EsinO, 



(1.2) 
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|<7<2)|sin0(2)^ 



ifi 2) ■ EsinO 

' ' = arcsm , 
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(4.15) 
(4.16) 

(4.17) 



such that the only free parameters are the incidence angle 6, the width of the barrier L and the strength of the 
interaction Dt- The transmission coefficient is then calculated numerically, using the definition Eq. p.7|) and the 
relation Eq. (|4.14|) . 

In order to investigate the Klein tunneling, including the effects of the K o K' mixing, in Figs. |3H5] we plot the 
transmission probability for a tensor barrier (blue line) against a 7*^ barrier given by Eq. p.9p (red line). For a lower 
value of the electron energy {E = 50 meV), we notice that for a relatively low barrier (half of the electron energy), 
the tensor barrier has a higher acceptance than the 7° barrier — this is mostly visible for a width of i = 110 nm, 
where the transmission coefficient is close to one for 6 < 38° for tensor and 9 < 17° for 7° barrier, as can be seen in 
Fig. El Increasing the height, the 7" barrier picks up more "structure", becoming transparent for central incidence 
and at several special angles, whereas the tensor interaction shows a less rich "structure" . For even higher barrier 
the behavior follows the same trend except that there is a smaller number of angles for which perfect tunneling is 
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FIG. 4: (Color online) Same as in Fig. O for an energy of the incoming electron E = 80 meV. The widths of the barrier are 
L = 110 nm (upper panel) and L — 50 nm (lower panel), and Dt denotes the height of the barrier. 
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FIG. 5: (Color online) Same as in Fig. [S] for an energy of the incoming electron E — 110 meV. The widths of the barrier are 
L = 110 nm (upper panel) and L = 50 nm (lower panel), and Dt denotes the height of the barrier. 



allowed. In the second set of plots (depicted in Fig. 2]) we increase the electron energy, while keeping the widths of the 
barriers constant. Choosing the electron energy identical to the value used in the calculation presented ini^ {E = 80 
meV), we find that the transmission probabilities for tensor and 7° barriers are almost identical. Finally, inspecting 
the last set of plots (Fig. [5|), we find that for even higher electron energies {E = 110 meV), the roles of the barriers 
are reversed: for a lower height, the tensor barrier has a lower acceptance compared with the 7° barrier, whereas for 
higher barriers the tensor component exhibits more angles for which perfect tunneling takes place. Our results show 
that the effects of the K o K' mixing, combined with the splitting of the wave inside the barrier, are more prominent 
for higher energies of the incoming electron. We speculate that this is a result of the fact that the exponential factors 
of the wave function inside the barrier, Eq. (|4.4p . depend on the absolute value of the electron energy E (and the 
barrier height Dt), whereas the angles Eq. (|4.17p only depend on the ration E/Dt- 
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V. SUMMARY AND CONCLUSIONS 

111 this paper, we have used the four component spinor formahsni to investigate the transmission coefficients for 7^, 
7°, /, 75 and 7^7^ one-dimensional barriers in graphene. We have worked in the framework of the gauge theoretical 
model proposed in Refs. |39||40| . where the electron-hole and carbon dynamics are simulated by a scalar and gauge 
fields. The different operators in the Dirac equation can be related to various types of defects in graphene: the 7^ 
operator is related to the breaking of the sublattice symmetry, 7"'^7^ describes a topological defect of the graphene 
sheet, whereas 7*^ can be associated with a chemical potential, and / and 75 appear in effective gauge theories for 
graphene^SiiS. 

For the 7° barrier we reproduce the results of the previous calculations performed within the two component spinor 
formalism^. For the 7^, / and 75 barriers we found that the transmittance has a strong dependence on the barrier 
height and incidence angle of the electron. In particular, under certain conditions these interactions generate an 
opaque barrier which gives rise to charge confinement. From an experimental point of view, the most prominent 
example is the breaking of the sublattice symmetry associated with the 7^ operator, which can be implemented via 
covalent bonding of foreign atoms to carbon atoms. Hence, according to our calculation, chemical bonding can be used 
as an instrument for controlling charge clustering in graphene. The remaining one-dimensional barrier, corresponding 
to the tensor interaction 7^7^ = iSa, exhibits an unusual feature. Specifically, it appears that the wave function of the 
incoming electron is divided into two plane waves inside the barrier, propagating at different angles and with different 
momenta. This translates into a distinct behavior of the transmission probability as a function of the incidence 
angle of the electron: while the Klein tunneling is observed just like in the case of the 7" barrieri^^i^, for particular 
values of the energy the scattering patterns of the two barriers can differ significantly (sec also Figs. [3HS1). We hope 
that understanding the nature of the tunneling states will be of further use in studies of the transport properties of 
graphene-based devices. 
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